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ABSTRACT 



We examine the evolution of the IGM Lya optical depth distribution using 
the transmitted flux probability distribution function (PDF) in a sample of 63 
QSOs spanning absorption redshifts 1.7 < z < 5.8. The data are compared to 
two theoretical r distributions: a model distribution based on the density distri- 
bution of Miralda-Escude et al. (2000) (MHR00), and a lognormal distribution. 
We assume a uniform UV background and an isothermal IGM for the MHR00 
model, as has been done in previous works where transmitted flux statistics have 
been used to infer an abrupt change in the IGM at z ~ 6. Under these assump- 
tions, the MHR00 model produces poor fits to the observed flux PDFs at redshifts 
where the optical depth distribution is well sampled, unless large continuum cor- 
rections are applied. However, the lognormal r distribution fits the data at all 
redshifts with only minor continuum adjustments. We use a simple parametriza- 
tion for the evolution of the lognormal parameters to calculate the expected mean 
transmitted flux at z > 5.4. The lognormal r distribution predicts the observed 
Lya and Ly/3 effective optical depths at z > 5.7 while simultaneously fitting the 
mean transmitted flux down to z — 1.6. In contrast, the best-fitting power-law 
under-predicts the amount of absorption both at z > 5.7 and at z < 2.5. If the 
evolution of the lognormal distribution at z < 5 reflects a slowly-evolving density 
field, temperature, and UV background, then no sudden change in the IGM at 
z ~ 6 due to late reionization appears necessary. We have used the lognormal 
optical depth distribution without any assumption about the underlying density 
field. If the MHR00 density distribution is correct, then a non-uniform UV back- 
ground and/or IGM temperature may be required to produce the correct flux 
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PDF. We find that an inverse temperature-density relation greatly improves the 
PDF fits, but with a large scatter in the equation of state index. The lognormal 
r distribution therefore currently offers the best match to the observed flux PDF 
and the most reliable predictor for the transmitted flux at high redshift. 

Subject headings: cosmology: observations — cosmology: early universe — in- 
tergalactic medium — quasars: absorption lines 

1. Introduction 

The Lyct forest serves as our most fundamental probe of the evolution of the intergalactic 
medium (IGM). While numerous models have been proposed for the underlying density 
field (see Rauch 1998 for a review), the current consensus is a self-gravitating network of 
filamentary structures collapsing out of initially Gaussian density perturbations. Given a 
description of the IGM that relates density and transmitted flux, one can compute various 
cosmological parameters and examine the large-scale evolution of the Universe. 

Perhaps the most dramatic inferences drawn from the evolution of Lya transmitted 
flux is that the reionization of the IGM may have ended as late as z ~ 6.2 (Becker et al. 
2001; White et al. 2003; Fan et al. 2002, 2006). This conclusion is based not only on the 
appearance of complete Gunn-Peterson troughs in the spectra of QSOs at z > 6, but on the 
accelerated decline and increased variance in the mean transmitted flux at z > 5.7 (Fan et 
al. 2006). Late reionization is potentially at odds with the transmitted flux seen towards the 
highest-redshift known QSO, SDSS Jl 148+5251 (z Q so = 6.42, White et al. 2003, 2005; Oh 
& Furlanetto 2005). In addition, the fact that the observed number density of Lya-emitting 
galaxies does not evolve strongly from z ~ 5.7 to z ~ 6.5 implies that the IGM is already 
highly ionized at these redshifts (Hu et al. 2004; Hu & Cowie 2006; Malhotra & Rhoads 2004, 
2006; Stern et al. 2005). Galactic winds (Santos 2004) or locally ionized bubbles (Haiman 
& Cen 2005; Wyithe & Loeb 2005; Furlanetto et al. 2004, 2006) may allow Lya photons to 
escape even if the IGM is significantly neutral. Additional arguments may be made about 
the thermal history of the IGM (Theuns et al. 2002; Hui & Haiman 2003) or the apparent 
size of the transmission regions around z ~ 6 QSOs (Mesinger & Haiman 2004; Mesinger et 
al. 2004; Wyithe & Loeb 2004; Fan et al. 2006). However, the evolution of the Lya forest 
remains the strongest evidence for late reionization. 

Still, the significance of the disappearance of transmitted flux at z ~ 6 has been highly 
debated (Songaila & Cowie 2002; Songaila 2004; Lidz et al. 2006a). As Songaila & Cowie 
(2002) pointed out, the mean transmitted flux in an inhomogeneous IGM will depend strongly 
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on the underlying density distribution, or more precisely, on the optical depth distribution. 
At z ~ 6, any transmitted flux will arise from rare voids, which lie in the tail of the optical 
depth distribution. Using a sample of 19 QSOs at z > 5.7, Fan et al. (2006) showed that the 
evolution the mean transmitted flux at z ~ 6 diverges significantly from that expected for 
a commonly- used model of the IGM density (Miralda-Escude et al. 2000, referred to herein 
as MHR00). The question, then, is whether the MHR00 model describes the distribution of 
optical depths accurately enough to make reliable predictions at very high redshift. 

In this paper we examine two theoretical optical depth distributions and their predictions 
for the Lja transmitted flux probability distribution function (PDF) The first is based on 
the gas density distribution given by MHR00, which has been used to make claims of late 
reionization. Their density distribution is derived from simple arguments about the dynamics 
of the IGM (see §3.2) and matches the output of an earlier numerical simulation (Miralda- 
Escude et al. 1996). In order to compute optical depths, assumptions must be made about 
the ionizing background and the thermal state of the IGM. As other authors have done, 
we will primarily consider a uniform UV background and an isothermal IGM. In §5 we will 
briefly generalize to a non-isothermal equation of state. 

The second case we consider is a simple lognormal optical depth distribution. This 
choice can be motivated in at least two ways. Initially Gaussian density perturbations will 
give rise to a lognormal density field when the initial peculiar velocity field is also Gaussian 
(Coles & Jones 1991). Indeed, Bi et al. (1992) demonstrated that a lognormal density 
distribution can produce many properties of the observed Lja forest (see also Bi et al. 1995; 
Bi & Davidsen 1997). More generally, however, a lognormal distribution naturally arises as 
a result of the central limit theorem when a quantity is determined by several multiplicative 
factors. For optical depth, these factors are density, temperature, and ionization rate. Here 
we will consider the lognormal distribution to be a generic distribution with the desirable 
properties of being non-zero and having a potentially large variance. Our main conclusions 
will not depend on any assumptions about the underlying density field. 

The transmitted flux PDF has been used to constrain a variety of cosmological param- 
eters (e.g., Rauch et al. 1997; Gaztahaga & Croft 1999; McDonald et al. 2000; Choudhury 
et al. 2001; Desjacques & Nusser 2005; Lidz et al. 2006b), with many authors assuming an 
optical depth distribution similar to one we consider here. We will examine the distributions 
themselves and their evolution with redshift by attempting to fit the models to the observed 
flux PDFs from a large sample of Keck HIRES data spanning Lja absorption redshifts 
1.7 < z < 5.8. We introduce the data in §2. In §3 the optical depth distributions are derived 
and used to fit the observed flux PDFs. We find that the lognormal distribution provides a 
better fit to the data at all redshifts where the optical depth distributions are well sampled. 
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In §4 we perform a simple fit to evolution of the lognormal distribution and use it to predict 
the mean transmitted flux at z > 5.7. In §5 we modify the model distribution by applying 
a non-isothermal equation of state. Finally, our results are summarized in §6. 

2. The Data 

Observations were made using the HIRES spectrograph (Vogt et al. 1994) on Keck I 
between 1993 and 2006. Targets are listed in Table I. QSOs at zq S o < 4.8 were observed 
using the original HIRES CCD and were reduced using the MAKEE package written by Tom 
Barlow. QSOs at zqso > 4.8 were observed using the upgraded detector and reduced using 
a custom set of IDL routines as described in Becker et al. (2006). The IDL package is based 
on the optimal sky subtraction technique of Kelson (2003). For nearly all of our observations 
we used an 0'/86 slit, which gives a velocity resolution FWHM of Av = 6.7 km s _1 . 

We will return to the issue of continuum fitting in §3.3. For now we will describe 
our baseline fitting procedure for quasars at various redshifts. For objects at zqso < 4.8, 
individual exposures were typically bright enough that a continuum could be fit to individual 
orders. This was done by hand using a slowly varying spline fit. The orders were then 
normalized prior to combining. At higher redshifts, we performed a relative flux calibration 
of each exposure using standard stars. The individual exposures were then combined prior 
to continuum fitting. A spline fit was again used for QSOs at zqso < 5.4. However, since 
the transmission regions at z > 5 rarely, if ever, reach the continuum, the fits were of a very 
low order and intended only to emulate the general structure of continua observed in lower 
redshift QSOs (e.g., Telfer et al. 2002; Suzuki 2006). For zqso > 5.7 we used a power law fit 
to the continuum of the form f v oc z/ -0 5 . 

Determining a quasar continuum is generally a subjective process whose accuracy will 
depend strongly on how much of the continuum has been absorbed (see Lidz et al. 2004b 
for a discussion). At z ~ 3, much of the spectrum will still be unabsorbed and errors in 
the continuum fit will depend on signal-to-noise of the data and the personal bias of the 
individual performing the fit. For high-quality data, errors in the continuum at z ~ 3 
should be < 1%. This uncertainty will increase with redshift as more of the continuum gets 
absorbed. By z ~ 5.5, very few transmission regions remain and the continuum must be 
inferred from the slope of the spectrum redward of the Lya emission line. However, the 
spectral slope may have an unseen break near Lya. In addition, echelle data are notoriously 
difficult to accurately flux calibrate. We therefore expect our power-law continuum estimates 
at z ~ 6 to be off by as much as a factor of two. 
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3. Flux Probability Distribution Functions 

3.1. Observed PDFs 

Observed transmitted flux probability distribution functions (PDFs) were taken from 
spectra of the 63 quasars listed in Table 1. In order to avoid contamination from the prox- 
imity region and from O Vl/Ly/3 absorption, we limited our analysis to pixels 10000 km s" 1 
redward of the Lya emission line and 5000 km s" 1 redward of the O VI emission line. The 
offsets were made intentionally large to account for possible errors in the QSO redshifts. 
In order for each region to contain enough pixels to be statistically significant yet avoid 
strong redshift evolution within a sample, we divided the Lya forest in each sightline into 
two sections covering ~ 60 A rest wavelength. Regions containing damped Lya systems 
were discarded. We further exclude wavelengths covered by the telluric A and B bands. 
Other atmospheric absorption due to water vapor was typically weak compared to the Lya 
absorption at the same wavelength and so was ignored. Table 1 lists the redshift interval for 
each region of the Lya forest we examine. 

Metal lines can be a significant contaminant in the Lya forest, particularly at lower 
redshifts. We therefore removed as many lines as could be identified either by damped Lya 
absorption or from multiple metal lines at the same redshift. In addition to the doublets 
C IV, Si IV, and Mg II, we searched for coincidences of Si II, Si ill, C II, O I, Fe II, Al II, 
and Al III. For exceptionally strong systems we also masked weaker lines such as Cr II, 
Ti II, S II, and Zn II. Lines in the forest were masked according to the structure and extent 
of lines identified redward of Lya emission. Very strong line that could be identified only 
from their presence in the Lya forest (e.g., saturated C IV) were also masked. However, we 
did not mask weak lines found in the forest without counterparts redward of Lya emission. 
Doing so would preferentially discard pixels with low Lya optical depth (where the metal 
lines can be seen), introducing a potentially larger bias in the PDF than the one incurred 
by leaving the contaminated pixels in the sample. In any case, our primary concern is with 
strong metal lines that could mimic saturated Lya absorption. Weak metal lines are not 
expected to significantly alter the flux PDF. 

The observed transmitted flux PDF for each region was computed in normalized flux 
bins of 0.02. Errors were computed using bootstrap resampling (Press et al. 1992). Each 
region was divided into many short sections spanning 200 km s _1 , and 1000 replicates of 
each region were constructed by randomly drawing sections with replacement. For this work 
we have used only the diagonal elements of the error matrix. As noted by McDonald et al. 
(2000) and Desjacques & Nusser (2005), ignoring the off-diagonal elements when performing 
X 2 fitting can have a significant effect on the width of the x 2 distribution, but has only 
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a small effect on the values of the best-fit parameters. For comparison, we have repeated 
the analyses presented in this paper using purely Poisson errors and have obtained nearly 
identical results. 



3.2. Theoretical PDFs 

We will examine two possible distributions for Lya optical depths: one based on the 
gas density distribution given by MHROO, and the other a lognormal distribution. In this 
section we derive the expected flux PDF for each case. 



3.2.1. MHROO model 

The MHROO gas density distribution is derived analytically based on assuming that the 
density fluctuations are initially Gaussian, that the gas in voids is expanding at constant 
velocities, and that the densities are smoothed on the Jeans length of the photoionized gas. 
The resulting parametric form for the volume-weighted density distribution is 



pMHR00 (A)=Aexp 



(A- 2 / 3 - Coy 

2(25 /3) 2 



(1) 



where A = p/p is the gas overdensity and A, Co, So, and (3 are constants. We take So = 
7.61/(1 + z) and (3 from Table 1 of MHROO, which produces good fits to ACDM (tt m = 0.4) 
simulation of Miralda-Escude et al. (1996). We then set A and C such that the total area 
under P^ HR00 (A) and the mean overdensity are both equal to one. Parameters for redshifts 
other than those listed in MHROO are linearly interpolated. 

To convert from densities to optical depths, assumptions must be made about the ion- 
izing background radiation and the thermal state of the gas. The Lya optical depth of a 
uniform IGM would be 

7T6 2 

t u = f a \ a H~ l (z)n ni , (2) 

m e c 

where f a is the Lya oscillator strength, X a = 1216 A, and H(z) is the Hubble constant at 
redshift z (Gunn & Peterson 1965). In the case of photoionization equilibrium, the optical 
depth t(A) for an overdensity A can be expressed in terms of the H I ionization rate T, and 
the recombination coefficient a as (Weinberg et al. 1997) 

T(A) K hT(A,zW ' ( ) 
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where a depends on the temperature as ct(T) oc T~°- 7 for T ~ 10 4 K (Abel et al. 1997). 
The IGM temperature will generally depend on the density, which is typically expressed as 
a power-law equation of state, T(A) = T A 1-7 (e.g., Hui & Gnedin 1997). However, as 
other authors have done, we will assume a uniform UV background and an isothermal IGM 
(Songaila & Cowie 2002; Songaila 2004; Fan et al. 2002, 2006). Following Fan et al. (2002), 
we can then express the optical depth as a function of density, 



r(A) = to 
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|r-i 2 (*)J 



(4) 



where T_ 12 is the H I ionization rate in units of 10 12 s 1 . For comparison to other works 
(McDonald k Miralda-Escude 2001; Fan et al. 2002, 2006), we take r = 82, although the 
normalization depends on the choice of cosmology. Equations (1) and (4) can then be used 
to determine the expected distribution of optical depths, 
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G = r 
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Finally, we can convert to the expected distribution of normalized fluxes, F = e T , 



,MHR00 



(F) 



2(-lnF)0 3 + 1 V 2 F 



cxp 
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(5) 



(6) 



(7) 



for < F < 1, otherwise. The distribution of fluxes at a particular z is then fully specified 
by the ionization rate r_i2. 



3.2.2. Lognormal t distribution 



For the lognormal optical depth distribution, we make no assumptions about the un- 
derlying density field, temperature, or ionization rate. As discussed above, a lognormal 
distribution can be motivated either from arguments about the evolution of an initially 
Gaussian density field (Coles & Jones 1991; Bi et al. 1992) or by the central limit theorem. 
Here we consider it to be a generic model that may plausibly describe the distribution of 
optical depths. The lognormal distribution is described by two parameters, \i = (lnr), and 
a, which is the standard deviation of lnr, 



^Lognormal 
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This gives an expected distribution of transmitted fluxes, 



pLognormal 
F 



(F) 



cxp 



(ln(-lnF) -nf 
2^ 



(9) 
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for < F < 1, otherwise. There are obvious similarities between the MHROO and lognormal 
distributions, which should not be surprising if they are both expected to at least roughly 
describe the data. We will examine the differences between the two cases more closely in 
§4.1. 



3.3. Fitting the observed PDFs 

In order to match the observed flux PDF, we must account for various imperfections 
in the data. The most important of these is noise in the flux measurements, which will 
smooth out the PDF and create pixels with F < and F > 1. We incorporate this effect 
by convolving the ideal flux PDFs given by equations (7) and (9) with a smoothing kernel 
constructed separately for each flux bin. (Numerically, the smoothing is performed on bins 
much narrower than those used for the final PDFs). The kernel for a particular bin is 
a weighted sum of Gaussian kernels whose widths and weights are determined from the 
distribution of formal flux errors of pixels in that bin. The result is typically a kernel with a 
narrow core to account for pixels with low noise, and an extended tail for noisier pixels. This 
allows us to fit regions of the Lya forest where the data quality is highly inhomogeneous. 

Errors in the continuum level and the flux zero point will also affect the the observed 
PDF. A change in the continuum will cause the observed PDF to be stretched or compressed 
in proportion to the flux level. An error in the zero point, which may result either from 
imperfect sky subtraction or from spurious counts (i.e., cosmic rays) improperly handled by 
the spectrum extraction or combination routines, will also stretch or compress the observed 
PDF from the low-flux end. In fitting the PDFs we consider two cases: first, where we 
assume there are no errors in either the continuum or the zero point, and second, where 
the continuum level and zero point are treated as free parameters. We define the preferred 
continuum and zero point levels to be those which, if applied to the data, would allow 
the theoretical distributions to produce the best. However, when performing the fits, the 
adjustments are applied to the models and not to the data. The continuum and zero points 
are treated independently, such that a change in the zero point does not require a change in 
the continuum, and visa versa. We do not allow zero point corrections at z < 3, where few 
pixels have zero flux. This was found to have no significant impact on the other parameters. 

The results of the x 2 minimization fitting are summarized in Tables 2 and 3 for the 
MHROO and lognormal cases, respectively. The best-fitting PDFs are plotted in Figures 1 
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through 5. For each region, we show the observed PDF along with the best-fitting theoret- 
ical PDFs in the cases where no continuum or zero point corrections are made and where 
the continuum and zero point are allowed to vary. At z > 5, the MHROO and lognormal 
distribution provide very similar fits. This is not a surprise since, at these redshifts, we are 
sampling the low-optical depth tail of both distributions. The differences in the distributions 
increase at lower redshift. At 3 < z < 5, the best-fit MHROO distribution significantly under- 
predicts the number of pixels with very low optical depth unless a continuum correction is 
applied. In contrast, the best-fit lognormal distributions provide a reasonable fit to the data 
at all redshifts, with or without a change in the continuum. Both models under-predict the 
number of saturated pixels in some cases, although the discrepancy tends to be much larger 
for the MHROO distribution. 

In Figure 6 we compare the minimum reduced x 2 values for both models in the case 
where the continua and zero points are held fixed. At z < 2, there is a roughly even 
divide between regions that are better fit by the MHROO distribution and those that prefer 
the lognormal distribution. However, in most instances where the MHROO distribution is 
preferred, the fit is relatively poor (x 2 > 2). At z > 5, the fits are mostly comparable, as 
noted above. For 3 < z < 4, the lognormal distribution provides a reasonable fit and is 
strongly preferred over the MHROO model. 

The fits improve for both models when the continua and zero points are allowed to 
vary. Most of this improvement is the result of the continuum corrections. The effect is 
particularly large for the MHROO distribution, which implies that the MHROO model tends 
to require that a significant continuum correction be applied to the data in order to produce 
a good fit. In Figure 7 we plot the reduced x 2 values for these more general fits. As was the 
case without continuum and offset adjustments, the two distributions produce comparable 
fits at z > 5. However, at all lower redshifts, the lognormal distribution is preferred. 

As noted above, even at low redshift, where extended regions of the spectrum have very 
little absorption, the continuum fit may be in error due to a combination of noise and the 
personal bias of the individual applying the fit. However, at z < 4, the continuum error 
should be less than a few percent for reasonably high signal-to-noise data. In Figure 8 we 
plot the continuum correction preferred for both distributions as a function of redshift. The 
MHROO model requires the continuum to steadily increase with redshift over the continuum 
drawn by hand in order to account for the lack of pixels predicted to lie near the continuum 
(i.e., pixels with very low optical depth). In contrast, the lognormal r distribution naturally 
accommodates fluxes near the continuum and does not require a large continuum correction 
for z < 4.5. At z > 5.4, the preferred continuum adjustment has a large scatter for both 
theoretical distributions, since nearly all pixels have significant optical depth. 
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In Figure 9 we show examples of the best-fit continua overlaid on the corresponding 
regions of the Lja forest. While the shape of a QSO continuum can be somewhat ambiguous 
when convolved with the response function of the instrument, no undue effort has been made 
to fit the continua across every transmission peak. The lognormal distribution fits the data 
well when the continua are near their intuitive values, while the MHROO model requires the 
continua to be substantially higher. Fitting QSO continua is an inherently uncertain task. 
However, even when the continuum is allowed to vary, the lognormal r distribution produces 
a better fit than the MHROO model. 

In contrast to our results, Rauch et al. (1997) and McDonald et al. (2000) found good 
agreement between the observed flux PDFs from some of the same sightlines used here 
and the predictions from a numerical simulation with a density distribution similar to the 
MHROO model. The reason for this appears to lie in their treatment of the continuum. Both 
works apply a strong correction to their simulated spectra by placing the continuum at the 
maximum transmitted flux level for each pass through the simulation box (lO/i -1 Mpc, or 
~ 30 A at z = 4). This is a much higher-order correction than we consider here. In addition, 
McDonald et al. (2000) group all pixels with flux F > 1 into their bin at F — 1. This 
disguises the shape of the observed PDF for pixels with low optical depth, particularly at 
z = 3 — 4. By fitting pixels at all fluxes, we remain sensitive to the shape of the PDF near 
F ~ 1. Applying a low-order continuum correction is therefore not sufficient to obtain a 
good fit for the MHROO distribution. However, this works well in the lognormal case. Much 
of the discriminating power in the flux PDF occurs at very low optical depths. Therefore, 
unless more reliable continuum fits can be made, the success of the MHROO model in this 
regime is at best unclear. 



4. Redshift Evolution of Optical Depth 

4.1. Lognormal Parameters 

We have shown that a lognormal distribution of optical depths provides a good fit to 
the observed Lya transmitted flux PDF at all redshifts 1.7 < z < 5.8. In this section we 
examine the evolution of the lognormal distribution and use it to predict the evolution of the 
mean transmitted flux at z > 6. In Figure 10 we plot the lognormal parameters fi and o as a 
function of z. Both parameters evolve smoothly with redshift, as should be expected if they 
reflect a slowly-evolving density field, UV background, and temperature-density relation. 
The increase in jj, and decrease in a with z can both be understood primarily in terms of 
the evolution of a self-gravitating density field. At earlier times, the density contrast in the 
IGM will be lower. This will tend to produce a higher volume-weighted median r, which is 
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given by e M , as well as a smaller logarithmic dispersion in r, which is given by a. Since we 
do not have an a priori model for how the lognormal parameters should evolve, for this work 
we choose the simplest possible parametrization. Excluding points at z > 5.4, where the 
lognormal parameters depend on highly uncertain continuum levels, a linear fit in redshift 
gives 

H( z ) = (-9.35 ± 0.17) + (1.79 ± 0.04)(1 + z) , (10) 

a(z) = (4.19 ± 0.16) - (0.46 ± 0.03)(1 + z) . (11) 

These fits are plotted as dashed lines in Figure 10. 

We can compare the evolution of the MHR00 and lognormal r distributions and their 
predictions for the transmitted flux PDF. In Figure 11 we plot fiducial r and flux distributions 
for 2 < z < 6. Parameters for the lognormal distribution are calculated from equations (10) 
and (11). For the MHR00 model, values for T_ 12 are chosen to be consistent with the fitted 
values in Table 2. The vertical dotted lines indicate the range of optical depths that can be 
measured with good data. At z = 2 we are primarily sensitive to the high-r tail in both 
distributions. At higher redshifts, the peaks of the distributions shift towards higher values 
of r until we are sampling only the end of the low-r tail at z — 6. 

Differences in the shape of the transmitted flux PDF are largest at 3 < z < 5, where 
P t (t) is well-sampled. The fact that the lognormal r distribution is most strongly favored 
at these redshifts suggests that it is more likely to be useful in making predictions for the 
distribution of transmitted flux at z > 6. An important feature of the lognormal distribution 
is that it narrows with redshift more rapidly than the MHR00 distribution. It therefore 
predicts fewer pixels with measurable transmitted flux at z ~ 6 than does the MHR00 
model with a slowly evolving UV background. 



4.2. Mean transmitted flux 

We can use the redshift evolution of the lognormal distribution to predict the the evo- 
lution of transmitted flux at z > 6. The mean transmitted flux will be given by 

(F) = [ FP F {F)dF. (12) 
Jo 

It is conventional to express the mean flux in terms of an effective optical depth r c g = 
— In (F). For a distribution of optical depths, r c ff will be smaller than the true mean optical 
depth. We show measurements of r c g for Lja from Songaila (2004) and Fan et al. (2006) in 
Figure 12. The dashed line shows the best-fitting power-law to their data at z < 5.5 from 
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Fan et al. (2006). The deviation of the data from the power-law at z > 5.7 has been cited 
as the primary evidence for an abrupt change in the ionizing background at z ~ 6. We also 
show r^j as predicted by the evolution of the lognormal r distribution given by equations 
(10) and (11) as a solid line. We emphasize that the lognormal parameters were fit only 
to measurements at z < 5.4. Even so, calculated from the lognormal distribution both 
better fits the data at z < 5 and predicts the upturn in at z > 5.7. In Figure 13 we 
include the lower-redshift measurements of Kirkman et al. (2005). The power-law under- 
predicts the amount of Lja absorption at z < 2.5, while the lognormal distribution matches 
all observations at 1.6 < z < 6.2. 

Stronger constraints on the ionization state of the IGM can be set using Ly/3, which 
is a weaker transition than Lja by a factor of 6.2. In the lognormal case, this produces a 
distribution of Ly[3 optical depths with the same a as Lja but with fj,p — /i a — In 6.2. We 
can then compute the expected mean flux in the Lj/3 forest at redshift z by multiplying the 
mean transmission resulting from Ly/3 absorption at z by the mean transmission resulting 
from Lja absorption at zp = (l + z)\p/\ a — l. We show the r^ s measurements from Songaila 
(2004) and Fan et al. (2006) in Figure 14. These are computed directly from the transmitted 
flux and have not been corrected for foreground Lja absorption. The dashed line again 
shows the best-fit power-law to the points at z < 5.5 from Fan et al. (2006). The solid line 
shows the lognormal prediction. Here again, despite the fact that we have not used any Lj(3 
measurements to determine the optical depth distribution, r^ s predicted in the lognormal 
case is a better fit to the data at z < 5 and follows the upturn in r^ s at z > 5.7. 

Our purpose here is not to fully characterize the evolution of transmitted flux at all 
redshifts. We have simply identified a distribution of optical depths that describes the 
observed distribution of transmitted fluxes better than the commonly used model. The fact 
that this distribution evolves smoothly with redshift, and that the same evolution describes 
changes in the Lya forest as well at z ~ 6 as it does at z ~ 3 strongly suggests that the 
disappearance of transmitted flux at z > 6 is due to a smooth evolution of IGM properties. 
The lognormal prediction for falls slightly below some of the lower limits of Fan et al. 
(2006) at z ~ 6, but the prediction does not take into account the expected scatter in the 
mean flux or any small deviation from our adopted linear redshift evolution of the lognormal 
parameters. The important point is that the evolution of the mean transmitted flux can be 
well described by a smooth evolution in the underlying optical depths. When sampling only 
the tail of the r distribution, as at z ~ 6, a slight change in the optical depths will produce 
a large change in the transmitted flux. 
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4.3. UV background 

Liu et al. (2006) recently demonstrated that a semi- analytic model based on a lognormal 
density distribution can reproduce the observed rise in r c g at z > 5.7. However, they invoke 
a UV background that declines rapidly with redshift, decreasing by a factor of ~ 11 from 
z = 3 to 5, and by a factor of ~ 7 from z = 5 to 6. We have not assumed that the lognormal 
r distribution used here arises directly from a lognormal density distribution. However, if 
we assume a uniform UV background and an isothermal IGM, than we can calculate the H I 
ionization rate by inverting equation (4) and averaging over all densities. Doing so gives 

r_ u = 0.05(1^)"^, (is) 

where (r 1 / 2 ) 2 = e ^+°" 2 / 4 5 and we have used the fact that (A) = 1. 

In Figure 15 we show T_i2 calculated for each fitted region along with the mean values 
in bins of redshift. For comparison, the best-fit values of T_i2 for the model distribution are 
also shown. The lognormal values are somewhat higher than the model values, which are 
in turn roughly consistent with previous measurements (McDonald & Miralda-Escude 2001; 
Fan et al. 2006). However, we do not require the strong evolution in T_ 12 given by Liu et al. 
(2006) for the lognormal model. Transforming from densities to optical depths depends on 
a number of factors, and we do not presume that the assumptions implicit in equation (13) 
are valid. We merely point out that a lognormal r distribution is consistent with a slowly 
evolving UV background. 

We can also calculate the mean volume-weighted neutral fraction, 

/ H : = (5.5 x 10- 5 ) Q^j (i^) 1 (1 + z)^ 2 (r) , (14) 

where we have used H{z) w # ^m /2 (l + zf/ 2 . The mean optical depth for the lognormal 
distribution will be (r) = e M+cr I 2 . Calculating \x and a from equations (10) and (11), this 
gives /hi = [1.0,1.2,1.9,4.0,11,20] x 10~ 5 for z = [2,3,4,5,6,6.5]. The mean optical depth 
will depend strongly on the high-r tail of the distribution, which is poorly constrained at 
z > 4. However, the disappearance of transmitted flux at z > 6 is at least consistent with a 
highly-ionized IGM. 

5. An inverse temperature-density relation? 

We have shown that the simplest transformation of the MHR00 gas density distribution 
to optical depths provides at best an uncertain fit to the observed distribution of transmitted 
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fluxes. However, there are several ways to modify the expected r distribution. Here we 
consider a non-isothermal temperature-density relation. From equation (3) we have r oc 
T _0 7 r^ 1 A 2 . We will address the general case where either T or T may depend on A. For a 
power law T°- 7 T oc A a , this gives 



where T_ 12 is now the H I ionization rate at the mean density, and the temperature at the 
mean density is included in r . For a uniform UV background, the equation of state index 
will be 7 = 1 + 1.43a. 

Not surprisingly, adding a degree of freedom significantly improves the fits for many of 
our Lya forest regions. The fitting results are summarized in Table 4, and a sample of the 
fits are shown in Figure 16. There is a large scatter in the best-fit a at all redshifts when the 
continuum and zero point are allowed to vary. However, the mean value (a) = —0.36 ± 0.45 
(sample variance) suggests that T 07 T increases towards lower densities. For a uniform UV 
background, this implies an equation of state T(A) oc A 7-1 with 7 m 0.5. An index < 1 
disagrees with previous measurements using the flux PDF (Choudhury et al. 2001; Lidz et al. 
2006b; Desjacques & Nusser 2005). However, those works typically considered only 7 > 1, 
which is expected following reionization if overdense regions experience more photoionization 
heating and less adiabatic cooling than underdense regions. Radiative transfer effects may 
create a complex temperature-density relation if underdense regions are reionized by a harder 
UV background than the dense regions near ionizing sources (Bolton et al. 2004). For the 
flux PDF, 7 < 1 allows for a lower T (typically by ~ 20%, see Figure 15), creating more 
saturated pixels, while at the same time maintaining a low r in low density regions. The 
necessary continuum corrections also decrease, although they are still roughly half of those 
needed in the case of a — 0. Of course, it is possible that we are not measuring the real 
equation of state, and that the added degree of freedom simply compensates for some other 
aspect of the model distribution. A more careful treatment of this problem will be reserved 
for future work. 



We have analyzed the Lya transmitted flux probability distribution in a high-resolution 
sample of 63 QSOs spanning the absorption redshift range 1.7 < z < 5.8. Our main goal has 
been to assess how well the theoretical optical depth distribution commonly used to measure 
the H I ionization rate describes the observed flux PDF. We find that the MHR00 model, 
under the assumptions of a uniform UV background and an isothermal IGM, produces a 




(15) 



6. 



Conclusions 



- 15 - 



poor fit to the observed flux PDF at all redshifts where the optical depth distribution is well 
sampled. This discrepancy eases only if large continuum corrections are applied. 

In contrast, a lognormal distribution of optical depths fits the data well with only minor 
continuum adjustments. The parameters of the lognormal distribution evolve smoothly with 
redshift, as expected for a slowly evolving IGM, and reflect both an increase in the mean r 
and a decrease in the relative scatter in r with redshift. We have performed simple linear 
fits to the lognormal parameters at z < 5.4. The mean transmitted flux calculated from 
these fits matches the observations at 1.6 < z < 5.7 better than the best-fitting power law 
(Fan et al. 2006). In addition, extrapolating the lognormal evolution to z > 6 predicts the 
observed upturn in both Lya and Ly/5 effective optical depths. This strongly suggests that 
if a slowly evolving density field, ionizing background, and IGM temperature are responsible 
for the evolution of the Lja forest at z < 5, then there is no reason to suspect a sudden 
change in the IGM at z ~ 6. 

We emphasize that we have used the lognormal distribution as a phenomenological 
description of the optical depths only, and that the distribution may not hold for optical 
depths that are outside the dynamic range of the transmitted flux. Other factors, such 
a non-isothermal IGM or variations in the UV background are likely to be important in 
deriving the optical depth distribution from the underlying density field. We have explored 
the possibility of a non- isothermal IGM in the context of the MHR00 model. The best 
fits tend to favor an inverse temperature-density relation, where temperature increases with 
density. This is contrary to typical expectations for the balance between photoionization 
heating and adiabatic cooling (Hui & Gnedin 1997), and may be an artifact of some other 
feature that causes the MHR00 model to disagree with the data. However, as Bolton et al. 
(2004) point out, radiative transfer effects may create a complex thermodynamic state in 
the IGM. If gas at a given density can have a range of temperatures and/or ionization rates, 
then a MHROO-like density distribution may give rise to a r distribution that is closer to 
lognormal. 

The largest source of uncertainty in fitting the flux PDFs remains the continuum level. 
Much of the disagreement between the MHR00 model and observed PDFs stems from the 
lack of pixels predicted to have very low optical depths at z > 3. This can be at least partially 
remedied by adjusting the continuum (see also McDonald et al. 2000). However, a dramatic 
change in the IGM would still be required to explain the observed lack of transmitted flux at 
z ~ 6.2 (e.g., Fan et al. 2002, 2006). Future observation of z > 4 gamma-ray bursts, whose 
continuum is a simple power law, may help to establish the correct flux PDF. For now, we 
have identified an optical depth distribution that both fits the data down to z — 1.6 and 
captures the evolution of the mean transmitted flux at z > 5.7. If the lognormal distribution 
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truly reflects aspects of the real optical depth distribution, then the motivation for late 
reionization may be greatly diminished. 

The authors would like to thank Martin Haehnelt for stimulating conversations, as 
well as Tom Barlow and Rob Simcoe for reducing much of the data. We especially thank 
the Hawaiian people for the opportunity to observe from from Mauna Kea. Without their 
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Table 1. Fitted Lya Forest Regions 



QSO 


2QSO 


(^abs) 


z abs 1 


2 abs X 


median flux 


SDKS 11 148-1-5251 


6.42 


5.614 


5.430 


5.802 


0.05 


snss iifno+0524 


6.30 


5.514 


5.339 


5.692 


0.11 


SDSS 11623-1-3112 


6.25 


5.522 


5.339 


5.709 


0.15 


SOSS 11048-1-4637 


6.23 


5.516 


5.339 


5.696 


0.15 


SDSS J0818+1722 


6.00 


5.590 


5.417 


5.766 


0.10 






5.221 


5.066 


5.417 


0.07 


SDSS J0002+2550 


5.82 


5.465 


5.339 


5.592 


0.12 






5.076 


4.910 


5.245 


0.11 


SDSS J0836+0054 


5.80 


5.455 


5.339 


5.573 


0.05 






5.067 


4.893 


5.245 


0.04 


SDSS J0231— 0728 


5.42 


5.043 


4.885 


5.206 


0.11 






4.730 


4.563 


4.885 


0.09 


SDSS J0915+4244 


5.20 


4.849 


4.707 


4.993 


0.07 






4.509 


4.373 


4.647 


0.07 


SDSS J1204-0021 


5.09 


4.747 


4.582 


4.887 


0.09 






4.428 


4.277 


4.582 


0.09 


SDSS J2225— 0014 


4.87 


4.513 


4.381 


4.647 


0.10 






4.234 


4.087 


4.381 


0.11 


BRI1202-0725 


4.69 


4.074 


3.929 


4.214 


0.10 


BRI2237— 0607 


4.56 


4.254 


4.126 


4.377 


0.07 


Q0246+1750 


4.44 


4.123 


3.988 


4.260 


0.05 






3.851 


3.716 


3.988 


0.07 


Q1055+4611 


4.15 


3.846 


3.719 


3.975 


0.02 






3.591 


3.460 


3.719 


0.03 


Q0000— 263 


4.13 


3.833 


3.704 


3.961 


0.05 






3.574 


3.447 


3.704 


0.05 


Ol 645-1-5520 


4.10 


3.798 


3.672 


3.927 


0.01 






3.543 


3.417 


3.672 


0.02 


BRI0241— 0146 


4.08 


3.779 


3.652 


3.906 


0.05 






3.523 


3.398 


3.652 


0.06 


Q0827+5255 


3.91 


3.623 


3.503 


3.748 


0.01 






3.389 


3.265 


3.503 


0.02 


Q0055-2659 


3.65 


3.381 


3.266 


3.499 


0.05 






3.149 


3.033 


3.266 


0.06 


Q1422+2309A 


3.63 


3.358 


3.243 


3.475 


0.02 






3.126 


3.011 


3.243 


0.02 


Q0930+2858 


3.44 


2.955 


2.845 


3.067 


0.07 


Q0642+44 


3.40 


2.927 


2.818 


3.037 


0.08 


Q0956+1217 


3.31 


3.061 


2.954 


3.169 


0.04 






2.845 


2.738 


2.954 


0.05 


HS0741+4741 


3.23 


2.772 


2.664 


2.876 


0.03 


Q0636+6801 


3.18 


2.931 


2.827 


3.036 


0.02 






2.720 


2.618 


2.827 


0.02 


Ql 140+3508 


3.16 


2.916 


2.813 


3.021 


0.03 






2.708 


2.605 


2.813 


0.03 


HS1011+4315 


3.14 


2.766 


2.657 


2.869 


0.04 
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Table 1 — Continued 



QSO 




\^-abs/ 




min 
lbs 


? max 
abs 


median flux error 


Q0449-1326 


3.10 


2.860 


2. 


757 


2.962 


0.04 






2.654 


2 


.552 


2.757 


0.07 


Q0940-1050 


3.08 


2.844 


2 


.743 


2.947 


0.04 






2.639 


2. 


538 


2.743 


0.05 


HS1946+7658 


3.07 


2.627 


2. 


.524 


2.728 


0.03 


Q2231-0015 


3.02 


2.780 


2. 


.680 


2.881 


0.05 






2.578 


2 


479 


2.680 


0.07 


Q1107+487 


2.98 


2.546 


2 


446 


2.646 


0.04 


Q1437+3007 


2.98 


2.547 


2 


.448 


2.648 


0.05 


Q0216+0803 


2.98 


2.574 


2. 


.487 


2.658 


0.15 


Q1437+3007 


2.98 


2.746 


2. 


648 


2.846 


0.04 


Q0216+0803 


2.98 


2.748 


2. 


.659 


2.843 


0.12 


Q1244+3133 


2.97 


2.541 


2 


439 


2.638 


0.10 


Q1511+0907 


2.89 


2.658 


2 


562 


2.756 


0.06 






2.464 


2 


.368 


2.562 


0.08 


Q1132+2243 


2.88 


2.652 


2 


.556 


2.750 


0.06 






2.456 


2. 


361 


2.556 


0.09 


HS0119+1432 


2.87 


2.643 


2 


.547 


2.740 


0.03 






2.452 


2 


.353 


2.547 


0.04 


Q1549+1919 


2.84 


2.613 


2. 


517 


2.707 


0.01 






2.419 


2. 


.324 


2.517 


0.01 


Q0528-250 


2.81 


2.595 


2 


.492 


2.683 


0.05 






2.398 


2. 


302 


2.492 


0.07 


Q2344+1228 


2.79 


2.374 


2. 


.280 


2.470 


0.09 


HS1700+6416 


2.74 


2.525 


2. 


.432 


2.619 


0.01 






2.339 


2. 


244 


2.432 


0.02 


Q1442+2931 


2.66 


2.264 


2. 


.169 


2.352 


0.03 


Q1009+2956 


2.65 


2.436 


2 


.345 


2.527 


0.02 






2.252 


2 


.162 


2.343 


0.02 


Q1358+1134 


2.58 


2.370 


2. 


.282 


2.461 


0.15 


Q2343+1232 


2.58 


2.190 


2. 


.101 


2.281 


0.10 


Q2206-199N 


2.57 


2.356 


2 


.269 


2.447 


0.03 






2.188 


2. 


105 


2.269 


0.04 


Q1623+2653 


2.53 


2.323 


2 


.235 


2.411 


0.06 






2.146 


2 


.058 


2.235 


0.10 


Q0841+1256 


2.51 


2.127 


2. 


038 


2.214 


0.12 


Q0237-233 


2.24 


2.050 


1 


966 


2.128 


0.06 


Q1225+3145 


2.21 


2.016 


1 


.938 


2.098 


0.03 






1.857 


1 


777 


1.938 


0.04 


Q0421+019 


2.05 


1.870 


1 


795 


1.947 


0.08 


Q0119-0437 


1.98 


1.807 


1 


.733 


1.876 


0.14 


Q0058+0155 


1.96 


1.797 


1 


734 


1.859 


0.12 



a Mean absorption redshift. 



Table 2. Best-Fit MHROO Model Parameters (Isothermal) 



QSO (^abs) a ^Vbin b Continuum and zero point fixed Continuum and zero point allowed to vary 
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3 17 
. 1 I 


u 


9fi 


1 


137 


017 


7^ 
U. i o 


SDSS T1 204 — 0021 


4 


428 


79 


o 


36 


6.64 


o 


28 


1 


194 


0.005 


1.60 


RRT9937— 0607 

£J A1ZZO I UUU / 


4 


254 


84 


Q 


64 


4.94 


o 


50 


1 


094 


-0.005 


1.62 


SDSS T9995 — 0014 


4 


234 


86 


o 


28 


5.13 


o 


21 


1 


170 


-0.002 


1.29 


O02464-1 750 


4 


123 


65 


o 


40 


2.55 


o 


34 


1 


069 


0.009 


1.11 


RRT1 909 — 0795 


4 


074 


87 


o 


29 


10.08 


o 


23 


1 


182 


0.011 


2.21 


O0246+1 750 


3 


851 


78 


o 


55 


5.26 


o 


39 


1 


073 


-0.005 


2.11 


D1 0554-461 1 


3 


846 


62 


o 


24 


4.68 


o 


22 


1 


091 


0.002 


1.60 


O0000—963 

V^UUUU ZDO 


3 


833 


70 


o 


34 


5.50 


o 


28 


1 


097 


0.010 


1.56 


01 6454-5590 


3 


798 


60 


o 


33 


5.49 


o 


31 


1 


076 


0.002 


1.53 


RRT0941 —01 46 


3 


779 


67 





40 


6.43 


o 


29 


1 


076 


0.019 


1.86 


008974-^9^ 


Q 

o 


fi93 

UZrO 


55 


u 


3 1 

o± 


4 Q3 


u 


^9 


1 


059 


000 


l.UU 




3 


591 


63 


o 


51 


7.66 


o 


41 


1 


052 


0.004 


1.73 


Q0000-263 


3 


574 


73 





44 


5.32 





32 


1 


048 


0.002 


3.33 


Q1645+5520 


3 


543 


61 





33 


4.27 





28 


1 


055 


0.000 


1.17 


BRI0241-0146 


3 


523 


74 





32 


4.96 





26 


1 


064 


0.000 


1.92 


Q0827+5255 


3 


389 


61 





33 


5.88 





27 


1 


060 


0.004 


1.17 


Q0055-2659 


3 


381 


67 





69 


3.02 





51 


1 


030 


0.010 


1.11 


Q1422+2309A 


3 


358 


58 





56 


4.41 





42 


1 


037 


0.008 


1.16 


Q0055-2659 


3 


149 


73 





54 


4.75 





44 


1 


014 


0.013 


4.44 


Q1422+2309A 


3 


126 


57 





44 


2.99 





37 


1 


025 


0.004 


0.89 


Q0956+1217 


3 


061 


64 





37 


4.47 





29 


1 


035 


0.001 


1.50 


Q0930+2858 


2 


955 


74 





49 


1.24 





44 


1 


Oil 




1.02 


Q0636+6801 


2 


931 


56 





50 


3.02 





39 


1 


017 




1.30 


Q0642+44 


2 


927 


76 





38 


2.12 





29 


1 


030 




0.95 


Q1140+3508 


2 


916 


60 





58 


4.64 





39 


1 


023 




2.05 


Q0449-1326 


2 


860 


65 





38 


1.60 





30 


1 


024 




0.54 


Q0956+1217 


2 


845 


66 





62 


3.04 





45 


1 


018 




1.75 


Q0940-1050 


2 


844 


61 





39 


3.14 





30 


1 


019 




1.81 


Q2231-0015 


2 


780 


64 





37 


3.12 





30 


1 


021 




1.92 
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Table 2 — Continued 



QSO ( z abs) a ^bin b Continuum and zero point fixed Continuum and zero point allowed to vary 











p 


1 o c 
-12 




r 


p_ 


-12 


V_'U11 L . U JJUllllr 


Xr 


HS0741+4741 


2 


772 


60 





49 


4 


59 





35 


1.024 


1.45 


HS1011+4315 


2 


766 


63 





44 


2 


46 





32 


1.019 


1.13 


Q0216+0803 


2 


748 


86 





21 


1 


30 





19 


1.010 


1.23 


Q1437+3007 


2 


746 


65 





49 


1 


14 





45 


1.005 


1.03 


Q0636+6801 


2 


720 


56 





35 


2 


75 





27 


1.017 


1.22 


Q1140+3508 


2 


708 


61 





56 


3 


12 





40 


1.014 


1.66 


Q1511+0907 


2 


658 


70 





42 


1 


74 





35 


1.011 


1.34 


Q0449-1326 


2 


654 


69 





39 


1 


18 





45 


0.990 


0.91 


Q1132+2243 


2 


652 


71 





39 


2 


72 





27 


1.022 


1.55 


HS0119+1432 


2 


643 


58 





54 


1 


34 





48 


1.005 


1.15 


Q0940-1050 


2 


639 


66 





38 


3 


85 





24 


1.029 


1.80 


HS1946+7658 


2 


627 


59 





36 


2 


45 





28 


1.015 


1.10 


Q1549+1919 


2 


613 


53 





50 


3 


14 





43 


1.011 


1.17 


Q0528-250 


2 


595 


65 





25 


5 


44 





16 


1.022 


4.54 


Q2231-0015 


2 


578 


74 





29 


1 


56 





26 


1.007 


1.52 


Q0216+0803 


2 


574 


89 





18 


1 


32 





19 


0.991 


1.31 


Q1437+3007 


2 


547 


73 





45 


1 


44 





41 


1.005 


1.43 


Q1107+487 


2 


546 


72 





55 


1 


05 





49 


1.005 


0.96 


Q1244+3133 


2 


541 


81 





22 


2 


34 





20 


1.011 


2.12 


HS1700+6416 


2 


525 


54 





57 


2 


41 





44 


1.007 


1.57 


Q1511+0907 


2 


464 


76 





31 


1 


03 





30 


1.002 


1.02 


Q1132+2243 


2 


456 


74 





43 


1 


28 





48 


0.994 


1.21 


HS0119+1432 


2 


452 


60 





32 


1 


46 





29 


1.005 


1.36 


Q1009+2956 


2 


436 


56 





48 


2 


05 





43 


1.004 


1.83 


Q1549+1919 


2 


419 


55 





68 


3 


77 





48 


1.008 


2.60 


Q0528-250 


2 


398 


72 





30 


1 


10 





34 


0.993 


0.99 


Q2344+1228 


2 


374 


79 





29 


2 


15 





28 


1.002 


2.17 


Q1358+1134 


2 


370 


86 





10 


3 


11 





16 


0.942 


0.85 


Q2206-199N 


2 


356 


60 





37 


1 


09 





33 


1.004 


1.02 


HS1700+6416 


2 


339 


55 





41 


1 


78 





41 


1.000 


1.89 


Q1623+2653 


2 


323 


67 





41 


1 


69 





36 


1.004 


1.55 


Q1442+2931 


2 


264 


60 





53 





99 





62 


0.996 


0.77 


Q1009+2956 


2 


252 


59 





36 





67 





33 


1.003 


0.62 


Q2343+1232 


2 


190 


81 





28 


1 


14 





23 


1.010 


1.05 


Q2206-199N 


2 


188 


61 





33 





97 





31 


1.002 


0.99 


Q1623+2653 


2 


146 


77 





29 


1 


30 





35 


0.990 


1.07 


Q0841+1256 


2 


127 


81 





14 


2 


42 





22 


0.970 


1.11 


Q0237-233 


2 


050 


67 





13 


5 


61 





28 


0.967 


1.30 


Q1225+3145 


2 


016 


62 





32 


2 


00 





30 


1.002 


2.01 


Q0421+019 


1 


870 


70 





49 


5 


17 





95 


0.979 


2.97 


Q1225+3145 


1 


857 


63 





25 


1 


71 





24 


1.001 


1.70 


Q0119-0437 


1 


807 


83 





30 


4 


33 





74 


0.963 


1.93 


Q0058+0155 


1 


797 


79 





29 


3 


75 





64 


0.966 


1.96 



a Mcan absorption rcdshift. 

b Number of flux bins over which fit was performed. 
C H I ionization rate, in units of 10 — 1 s _1 . 
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d Factor by which to multiply the continuum in order for the model to produce the best fit. 
°Flux zero point that would allow the model to produce the best fit. 



Table 3. Best-Fit Lognormal Parameters 



QSO (^abs) a ^Vbin b Continuum and zero point fixed Continuum and zero point allowed to vary 













(j 


\ 


'2 




<7 


Cont.° 


Zero pt. 1 


X 


2 

r 


SDSS J1148+5251 


5 


614 


49 


1.81 





86 


1 


10 


2.21 


1 


24 





774 


0.009 


o 


49 


SDSS T0818+1722 


5 


590 


83 


2.71 


1 


83 


1 


98 


2.59 


1 


63 


1 


161 


0.002 


1 


91 


SDSS J1623+3112 


5 


522 


76 


1.58 





80 





83 


1.56 





75 


1 


088 


-0.001 


o 


87 


SDSS Tl 048+4637 


5 


516 


88 


1.58 


1 


19 


1 


42 


1.40 





83 


1 


427 


-0.010 


o 


99 


SDSS J1030+0524 


5 


514 


83 


1.53 


1 


07 


o 


97 


1.63 


1 


16 





983 


0.006 


o 


92 


SDSS J0002+2550 


5 


465 


68 


1.61 





79 


1 


21 


1.82 





94 





940 


0.009 


1 


14 


SDSS T0836+0054 


5 


455 


55 


1.54 


1 


10 


o 


64 


1.66 


1 


36 





866 


0.004 


o 


17 


SDSS J0818+1722 


5 


221 


73 


1.19 


1 


15 


o 


88 


1.29 


1 


31 





937 


0.007 


o 


71 


SDSS J0002+2550 


5 


076 


75 


1.19 





95 


1 


90 


1.42 


1 


25 





884 


0.019 


1 


43 


SDSS J0836+0054 


5 


067 


56 


1.29 


1 


00 


o 


72 


1.29 


1 


10 





922 


0.000 


o 


59 


SDSS J0231— 0728 


5 


043 


78 


0.80 


1 


23 


o 


91 


0.89 


1 


35 





974 


0.012 


o 


72 


SDSS J0915+4244 


4 


849 


71 


1.12 


1 


49 


1 


13 


1.28 


1 


61 


1 


004 


0.011 


o 


48 


SDSS J1204— 0021 


1 


747 


75 


0.26 


1 


36 


4 


23 


0.58 


1 


49 


1 


052 


0.033 


o 


99 


SDSS T0231 —0728 


4 


730 


77 


0.63 


1 


17 


1 


22 


0.64 


1 


38 





916 


0.009 


o 


96 


SDSS J2225— 0014 


4 


513 


80 


0.40 


1 


65 


o 


95 


0.49 


1 


63 


1 


028 


0.012 


o 


63 


SDSS J0915+4244 


4 


509 


77 


0.26 


1 


40 


1 


60 


0.40 


1 


43 


1 


036 


0.014 


o 


84 


SDSS J1204-0021 


4 


428 


79 


0.23 


1 


91 


1 


28 


0.31 


1 


87 


1 


026 


0.009 


o 


98 


RRT2237— 0607 


4 


254 


84 


-0.60 


1 


98 


1 


67 


-0.60 


1 


92 


1 


004 


0.008 


1 


64 


SDSS J2225— 0014 


4 


234 


86 


0.34 


1 


82 


o 


71 


0.35 


1 


88 





991 


0.004 


o 


69 


Q0246+1750 


4 


123 


65 


-0.35 


1 


56 


1 


54 


-0.32 


1 


51 


1 


016 


0.002 


1 


49 


BRI1202— 0725 


1 


074 


87 


-0.08 


2 


04 


2 


19 


0.05 


1 


98 


1 


032 


0.018 


1 


24 


Q0246+1750 


3 


851 


78 


-0.74 


2 


16 


1 


36 


-0.79 


2 


26 





987 


0.000 


1 


26 


01055+461 1 


3 


846 


62 


-0.15 


2 


02 


o 


89 


-0.12 


1 


96 


1 


019 


0.003 


o 


58 


Q0000— 263 


3 


833 


70 


-0.59 


1 


93 


1 


94 


-0.45 


1 


94 


1 


018 


0.010 


1 


19 


Q1645+5520 


3 


798 


60 


-0.54 


2 


06 


1 


15 


-0.53 


1 


96 


1 


022 


0.002 


o 


76 


BRI0241-0146 


3 


779 


67 


-0.78 


1 


85 


2 


39 


-0.61 


1 


99 


1 


004 


0.020 


1 


07 


Q0827+5255 


3 


623 


55 


-0.58 


2 


30 





80 


-0.65 


2 


20 


1 


015 


0.000 


o 


60 


Q1055+4611 


3 


591 


63 


-1.31 


2 


15 


1 


23 


-1.24 


2 


08 


1 


008 


0.003 


o 


96 


Q0000-263 


3 


574 


73 


-0.75 


2 


25 


1 


57 


-0.82 


2 


53 





981 


0.008 


1 


05 


Q1645+5520 


3 


543 


61 


-0.85 


2 


07 





76 


-0.85 


1 


98 


1 


009 


0.002 





70 


BRI0241-0146 


3 


523 


74 


-0.82 


2 


07 


1 


76 


-0.83 


2 


10 





997 


0.001 


1 


80 


Q0827+5255 


3 


389 


61 


-1.22 


1 


97 


1 


58 


-1.13 


1 


87 


1 


018 


0.004 





88 


Q0055-2659 


3 


381 


67 


-1.88 


2 


04 


1 


17 


-1.92 


2 


18 





994 


0.008 


1 


05 


Q1422+2309A 


3 


358 


58 


-1.79 


2 


03 


1 


56 


-1.62 


2 


16 


1 


005 


0.007 





82 


Q0055-2659 


3 


149 


73 


-1.67 


2 


48 


3 


41 


-2.16 


3 


32 





970 


0.011 


1 


15 


Q1422+2309A 


3 


126 


57 


-1.95 


1 


98 


1 


05 


-1.89 


1 


87 


1 


005 


0.002 





97 


Q0956+1217 


3 


061 


64 


-1.85 


2 


16 


1 


63 


-1.83 


2 


07 


1 


004 


0.002 


1 


67 


Q0930+2858 


2 


955 


74 


-2.28 


1 


99 


1 


05 


-2.44 


2 


15 





988 







89 


Q0636+6801 


2 


931 


56 


-2.23 


2 


13 





78 


-2.25 


2 


17 





999 







79 


Q0642+44 


2 


927 


76 


-2.00 


2 


20 





72 


-2.10 


2 


31 





992 







67 


Q1140+3508 


2 


916 


60 


-2.32 


2 


40 





89 


-2.35 


2 


43 





998 







93 


Q0449-1326 


2 


860 


65 


-2.16 


2 


07 





80 


-2.16 


2 


06 


1 


001 







83 


Q0956+1217 


2 


845 


66 


-2.65 


2 


54 


1 


15 


-2.86 


2 


81 





991 







78 


Q0940-1050 


2 


844 


61 


-2.07 


2 


28 


1 


31 


-2.22 


2 


54 





990 







93 


Q2231-0015 


2 


780 


64 


-2.25 


2 


03 


1 


48 


-2.20 


2 


00 


1 


003 




1 


50 



-25 - 



Table 3 — Continued 



QSO 


<^abs} a 


N hin h 


Continuum and 


zero point fixed 


Continuum and zero point allowed 


;o vary 










U 


xi 




U 


Cont. e Zsro pt.^ 


xi 


TTQ07/I 1 1 A 1 A 1 

xloU f 41+4 1 41 


O 770 

z. / / z 


OU 


-2.57 


2.36 


11/1 
1.14 


-2.46 


2.22 


1.006 


O (17 

u.y / 


liolUll+4olO 


Z. /DO 


KQ 
DO 


-2.39 


2.25 


n ki 
U.Dl 


-2.45 


2.33 


0.997 


o re 
U.Oo 


^JUZlDH-UoUo 


O 1 A O 
Z. ( 45 


OO 


-1.65 


1.91 


1 oo 
l.ZZ 


-1.99 


2.22 


0.971 


O OK 

u.yo 


ni /iQ7 i oon7 
^J14o ( +oUU f 


O 7/1 £ 
Z. / 4D 


K K 
DO 


-2.65 


1.96 


1.(0 


-2.89 


2.29 


0.988 


1 O/l 
1.Z4 


^UDoD-hDoUl 


o 7on 

Z. I ZU 


OD 


-2.26 


2.11 


n on 
u.yu 


-2.28 


2.14 


0.999 


O Ol 

u.y i 


^ 1 14U-)-oOUo 


O 7Hfi 
Z. i Uo 


K1 
Dl 


-2.65 


2.41 


n i a 
U. / 4 


-2.81 


2.60 


0.994 


n c;o 
u.oz 


m en i i non7 

^loii+uyu r 


O ft tTQ 
Z.DOo 


7(1 
IK) 


-2.62 


2.21 


1 /in 
1.4U 


-2.93 


2.58 


0.987 


o (if; 

u.yo 


^U44y — loZD 


O fi^/l 
Z.004 


KO 

Dy 


-2.52 


1.60 


o on 
z.zu 


-3.13 


2.32 


0.974 


i m 

1.1U 


HI 1 QO 1 OO/I Q 

^lloZ-|-ZZ4o 


o a eo 


/ 1 


-2.51 


2.43 


n oo 

u.yz 


-2.66 


2.59 


0.993 


o on 
U.oU 


rioUliy-rl4oZ 


O K/1 Q 

Z.D4o 


eo 
Oo 


-2.82 


2.08 


1 Q(l 

l.oU 


-3.13 


2.40 


0.991 


O 7Q 
U. ( O 


nnnyin 1 ne;n 
^JUy4U— lUoU 


O (JQO 

z.ooy 


kk 
DO 


-2.35 


2.62 


n kk 
U.DD 


-2.43 


2.73 


0.996 


O KQ 
U.DO 


TTQIO/lft I 7fiKQ 
rloiy4D+ f ODo 


O C07 
Z.DZ ( 


en 

oy 


-2.46 


2.12 


n 7(i 


-2.51 


2.17 


0.998 


O QO 

U.oU 


L<jio4y-i-iyiy 


O K1 Q 
Z.Dlo 


Oo 


-3.02 


2.14 


1 in 
1.1U 


-2.92 


2.01 


1.003 


O 07 

u.y * 


nriKOQ orn 

^JUOZo — ZOU 


O c:(ic; 

z.oyo 


K c; 
DO 


-1.68 


2.41 


O QC 

Z.oO 


-2.11 


2.83 


0.982 


1 KA 

l.OU 


nooQi nm e. 

t^ZZol — UU10 


o c;7o 

Z.O ( o 


7/1 
( 4 


-2.37 


2.26 


o nn 

Z.UU 


-2.87 


2.85 


0.975 


O 07 

u.y / 


l^JUZ10-|-UoUo 


O e.7/1 
Z.O f 4 


fiO 

oy 


-1.90 


1.97 


o no 
z.uz 


-2.83 


2.98 


0.939 


o flt; 

u.yo 


m a Q7 i qoo7 


O c:/l7 
Z.04 ( 


79 
( O 


-3.00 


2.34 


1 O/l 

l.y4 


-3.49 


3.01 


0.982 


o oo 

u.yo 


H1 1 CiT i /1C7 


O z^AR 
Z.04D 


70 

( z 


-3.17 


2.27 


117 
1.1 ( 


-3.55 


2.61 


0.989 


O 70 

u. / z 


H1 O/l /I 1 Q1 QQ 

^JlZ44H-oloo 


O e; a 1 
Z.041 


Q1 

ol 


-2.16 


2.46 


O 1 Q 

Z. lo 


-2.78 


3.15 


0.966 


O 7/1 

U. / 4 


liol (UU+041D 


O c:oc; 
Z.OZO 


C/1 
04 


-3.08 


2.38 


1 1 Q 

l.lo 


-3.31 


2.81 


0.995 


O ft7 
U.D ( 


m ei i i nnnv 

1^1011-t-uyu i 
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a Mcan absorption redshift. 

b Number of flux bins over which fit was performed. 
c Lognormal parameter fi = (lnr). 
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d Lognormal parameter a = stddev (lnr). 

°Factor by which to multiply the continuum in order for the model to produce the best fit. 
f Flux zero point that would allow the model to produce the best fit. 



-27- 



Table 4. Best-Fit MHROO Model Parameters (Non-Isothermal) 
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Table 4- 


-Continued 
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( z a,bs) s ' N bin b Continuum and zero 


point fixed 


Continuum and zero point allowed to vary 
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a Mcan absorption redshift. 

b Number of flux bins over which fit was performed. 
C H I ionization rate, in units of 10~ 1 s — 1 . 
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d Power-law index for the generalized temperature-density relation T°- 7 T oc A". 

°Factor by which to multiply the continuum in order for the model to produce the best fit. 

f Flux zero point that would allow the model to produce the best fit. 
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Fig. 1. — Fits to the Lya flux probability distribution functions for QSOs in our sample. Each 
set of panels is labeled with the QSO name and the mean absorption redshift. Histograms 
show the observed PDF. For each region, MHROO model fits assuming an isothermal IGM 
are shown on the left-hand side (red lines), while fits based on a lognormal r distribution are 
shown on the right-hand side (blue lines). Dotted lines indicate the best fit without adjusting 
either the continuum or the zero point. Solid lines show the best fits when the continuum 
and zero point are allowed to vary. The lognormal r distribution generally produces a good 
fit without large adjustments to the continuum, whereas at 3 < z < 5 the MHROO model 
fits tend to be poor unless a significant continuum adjustment is made. 
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Fig. 2. — Fits to the Lya flux probability distribution functions for QSOs in our sample, 
continued from Figure 1. 




Fig. 3. — Fits to the Lya flux probability distribution functions for QSOs in our sample, 
continued from Figure 1. 




Fig. 4. — Fits to the Lya flux probability distribution functions for QSOs in our sample, 
continued from Figure 1. 




Fig. 5. — Fits to the Lya flux probability distribution functions for QSOs in our sample, 
continued from Figure 1. 
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Fig. 6. — A comparison of the reduced \ 2 values for the best-fitting MHROO and lognormal 
t PDFs when the continuum and zero point are held fixed. Symbols indicate the mean 
absorption redshift of the fitted region of the Lya forest. The dotted line indicates where 
xl mhroo = X 2 r , Lognormal- Roughly half of the Lya regions at z < 3 are better fit by the 
MHROO PDF. Otherwise, the lognormal PDF is preferred. 
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Fig. 7. — A comparison of the reduced \ 2 values for the best-fitting MHROO and lognormal 
t PDFs when the continuum and zero point are allowed to vary. Note the change in scale 
from Figure 6. Symbols indicate the mean absorption redshift of the fitted region of the 
Lya forest. The dotted line indicates where Xr,MHRoo = Xr, Lognormal- The lognormal PDF is 
preferred at all redshifts, particularly at z < 5. 
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Fig. 8. — Continuum adjustments to the data that are required so that the observed distri- 
bution of Lya fluxes are best fit by the theoretical distributions. Squares show the contin- 
uum adjustment needed for the MHROO model. Diamonds show the continuum adjustment 
needed for the lognormal r distribution. The MHROO model value for SDSS J0818+1722 
at (^abs) = 5.590 lies outside the plot range, as indicated by the arrow. At z < 4, only 
the values for regions with median flux error < 0.05 are shown. The MHROO distribution 
requires a steadily increasing continuum adjustment with redshift to account for the lack of 
pixels predicted to lie near the continuum. At z > 5.4 the best-fitting continuum has a large 
scatter for both distributions due to how little transmitted flux remains in the forest. 
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Fig. 9. — Four examples of continuum adjustments needed so that the observed flux PDFs 
are best fit by the theoretical distributions. Each panel shows a sample of the Lya forest 
taken from the fitted region indicated by the QSO name and mean absorption redshift. The 
solid and dashed horizontal lines shows the continuum levels best fit by the MHROO and 
lognormal r distributions, respectively. The spectra have been binned to 13 km s _1 pixels 
for clarity. 
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Fig. 10. — Parameters for the lognormal distribution of Lya optical depths that produce the 
best fits to the observed flux PDFs. Here, \x = (hir) and a = stddev(lnr). The dashed 
lines show the best linear fits from equations (10) and (11). 
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Fig. 11. — Redshift evolution of the theoretical Lya optical depth and transmitted flux dis- 
tributions. The top panels show the r distributions for the indicated redshifts and ionization 
rates. Bottom panels show the corresponding transmitted flux PDFs. Distributions for the 
MHROO model are shown in red. Distributions for the lognormal r model are shown in 
blue. Parameters for the lognormal r distribution were calculated from fits to \i and a as a 
function of redshift (cf. equations 10 and 11). Vertical dotted lines indicate optical depths 
corresponding to 98% and 2% transmitted flux. The clearest differences in the predicted 
shapes of the flux PDFs occur at 3 < z < 5. The lognormal r distribution, which produces 
better fits to the data, narrows with redshift more rapidly than the MHROO r distribution. 
Hence, fewer pixels with measurable transmitted flux at z = 6 are predicted in the lognormal 
case. 
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Fig. 12. — The evolution of Lya effective optical depth with redshift, where r c g = — In (F). 
Data points are from Songaila (2004) (small circles) and Fan et al. (2006) (large circles and 
arrows, with colors matching their Figure 2). The dashed line shows the best-fit power-law to 
at z < 5.5 from Fan et al. (2006). The solid line shows calculated from the lognormal 
distribution of Lya optical depths, for which the parameters were fit at z < 5.4. A simple 
evolution in the lognormal r distribution predicts the upturn in at z > 5.5 and produces 
a better fit to the observed at 4 < z < 5. 
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Fig. 13. — Same as Figure 12, with on a logarithmic scale. We have also included lower- 
redshift measurements calculated from Kirkman et al. (2005), which exclude absorption from 
metal lines, Lyman limit systems, and damped Lya systems. The Kirkman et al. (2005) 
points are plotted as open squares with errors in the mean measurements. The power-law fit 
from Fan et al. (2006) (dashed line) under-predicts the amount of Lya absorption both at 
z > 5.7 and at z < 2.5. In contrast, calculated from the lognormal r distribution (solid 
line), provides a simultaneously good fit to all points at 1.6 < z < 6.2. 
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Fig. 14. — The evolution of Ly/3 effective optical depths with redshift. Data points at from 
Songaila (2004) (small circles) and Fan et al. (2006) (large circles and arrows, with colors 
matching their Figure 3). The data have not been corrected for foreground Lja absorption. 
The dashed line shows the best-fit power-law to at z < 5.5 from Fan et al. (2006). 
The solid line shows t^ s predicted purely from the lognormal distribution of Lya optical 
depths. Even though no independent fitting of Ly/3 fluxes was performed, the lognormal r 
distribution captures the upturn in at z > 5.5 and produces a better fit to the observed 
tL at 4 < z < 5. 
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Fig. 15. — H I ionization rate as a function of redshift. Small diamonds show T_i2 calcu- 
lated from the lognormal fits to individual regions, assuming a uniform UV background (cf. 
equation 13). Points fit to high-S/N data are shown in blue. Filled circles show the mean 
T_i2 for the lognormal model in redshift bins of 0.5, starting at z — 2. Vertical error bars 
show the standard deviation of points within a bin. Horizontal error bars show the range of 
redshift covered by all points within that bin. Squares show the mean T_i2 from the MHR00 
model fits for an isothermal IGM and uniform UV background (a = 0). Triangles show the 
mean T_ 12 from the MHR00 model fits when a is allowed to vary. 
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Fig. 16. — Examples of transmitted flux PDFs where the MHROO model fit is significantly 
improved by allowing a non-isothermal temperature-density relation, or more generally, 
T 0,7 r oc A a . Each set of panels is labeled with the QSO name and the mean absorp- 
tion redshift. Histograms show the observed PDF. For each section, MHROO model fits 
with a = are shown on the left-hand side (red lines), while fits with a treated as a free 
parameter are shown on the right-hand side (green lines). Dotted lines indicate the best 
fit without adjusting either the continuum or the zero point. Solid lines show the best fits 
when the continuum and zero point are allowed to vary. The mean value of a for all regions 
is (a) ~ —0.4. This may indicate that the UV background decreases with density, or that 
there exists an inverse temperature-density relation. Alternatively, finding a < may be an 
artifact of some other features of the MHROO model that causes it to disagree with the data. 



